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^ I The main result of this article is that we show that from supersymmetry we 

^1 can generate new superintegrable Hamiltonians. We consider a particular 

■ case with a third order integral and apply the Mielnik's construction in 
' supersymmetric quantum mechanics. We obtain a new superintegrable 

potential separable in Cartesian coordinates with a quadratic and quintic 

integrals and also one with a quadratic integral and an integral of order 

^ I seven. We also construct a superintegrable system written in terms of the 

\Q • fourth Painleve transcendent with a quadratic integral and an integral of 

'nI" . order seven. 
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'p. ; 1 Introduction 

> . 

^ . Superintegrability [1-14] and supersymmetric quantum mechanics 

^ ! (SUSYQM) [15-21] have attracted a lot of attention in recent years. Both of 

these fields have important applications in quantum chemistry, atomic 
physics, molecular physics, nuclear physics and condensed matter physics. 
Although they are two separate issues, many quantum systems such as the 
harmonic oscillator, the hydrogen atom and the Smorodinsky-Winternitz 
potential are both superintegrable and supersymmetric [21]. 
Superintegrability with third order integrals was the object of a series of 
articles [22-26]. The systems studied have a second and a third order 
integrals. They were studied by means of cubic and deformed oscillator 
algebras. The supersymmetric quantum mechanics approach was used [25] 
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and also higher order supersymmetric quantum mechanics [26] in order to 
calculate energies and wave functions. These articles indicate that 
superintegrability is closely connected to supersymmetry. We will show in 
this article that supersymmetry can provide a method of generating new 
superintegrable systems. We will consider two-dimensional systems 
separable in Cartesian coordinates. The separability implies the existence of 
a second order integral of motion. 

Let us recall some definitions concering superintegrability and 
supersymmetry. In classical mechanics a Hamiltonian system with 
Hamiltonian H and integrals of motion Xa 

H =^gikPiPk + V{x,p), Xa = fa{x,p), a = l,...,n-l , (1.1) 

is called completely integrable (or Liouville integrable) if it allows n 
integrals of motion (including the Hamiltonian) that are well defined 
functions on phase space, are in involution {H,Xa}p = 0, {Xa,Xi,}p = 0, 
a,b=l,...,n-l and are functionally independent ({, }p is a Poisson bracket). 
A system is superintegrable if it is integrable and allows further integrals of 
motion Yb{x,p), {H,Yb}p = 0, b=n,n+l,...,n+k that are also well defined 
functions on phase space and the integrals { if, Xi, ...,Yn^k} are 

functionally independent. A system is maximally superintegrable if the set 
contains 2n-l such integrals. The integrals Yf, are not required to be in 
evolution with Xi,...Xn-i, nor with each other. The same definitions apply 
in quantum mechanics but {H, Xa,Yh} are well defined quan- 
tum mechanical operators, assumed to form an algebraically independent set. 

In Section 2, we recall definitions and results of supersymmetric quantum 
mechanics. We also discuss some results obtained by B.Mielnik [27]. 
B.Mielnik showed that the factorization of second order operators is not 
necessarily unique. Supersymmetric quantum mechanics allows to find the 
eigenfunctions, the energy spectrum, creation and annihilation operators. In 
Section 3, we will consider a two-dimensional Hamiltonian consisting of two 
one-dimensional Hamiltonians that are superpartners. Such systems are by 
construction separable in Cartesian coordinates so a second order integral 
exists. From the creation and annihilation operators of the one dimensional 
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part we can generate a higher order integral of motion. The system is thus 
superintegrable. We show how these results allow us to recover known 
superintegrable systems with a third order integral that are special cases of 
a Hamiltonian written in terms of the fourth Painleve transcendent. In 
Section 4, we consider a particular case with a third order integral, apply 
the Mielnik's method and obtain a new superintegrable potential separable 
in Cartesian coordinates with a quadratic and quintic integrals and also one 
with a quadratic and seventh order integrals. We also constuct a 
superintegrable system written in terms of the fourth Painleve transcendent 
with a quadratic and seventh order integrals. 



2 Supersymmetry and Mielnik's factorization 
method 

We begin this Section by recalling definitions and results of supersymmetric 
quantum mechanics. We define two first order operators 

We consider the following two Hamiltonians which are called 
"superpartners" 

(2.2) 

There are two cases. The first is A^^^ 7^ 0, E^^^ ^ 0, A^ip^^^ ^ and 
^ 0. We have 

Ei^^ = Ei'^>0, ^f = ^A^i^), ^« = ^AVi^) . (2.3) 

and the two Hamiltonians are isospectral. This case corresponds to broken 
supersymmetry. 

For the second case the supersymmetry is unbroken and we have A-^q^^ = 0, 
^J^) = 0, A^ifj}^^ ^ and E^^^ ^ 0. Without lost of generality we take Hi as 
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having a zero energy ground state. We have 



(2.4) 

We can define the matrices 
They satisfy the relations 

[ij,Q] = [/j,Qt] = o, {g,Q} = {Q^gt} = o, {Q,gt} = // . (2.6) 

The operators Q, Q'^ are called "supercharges". We have a sl(l|l) 
superalgebra and Hi and H2 are superpartners. Supersymmetric quantum 
mechanics allow us to obtain the creation and annihilation operators. The 
operators M'^ and M with h'^ and h respectively the creation and 
annihilation operators for the Hamiltonian Hi 

M = A^bA, = A^b^A, (2.7) 

are thus the creation and annihilation operators for the Hamiltonian 7^2- 
Supersymmetric quantum mechanics with higher order supercharges has 
been studied [28-32]. The case with second order operators of the form 

= d'^ -2h{x)d + c{x), M = d"^ + 2h{x)d + c{x) . (2.8) 

was investigated. The case with a first and second order supersymmetry 
was also treated. 

As far as we could find, the generalized ladder operators appeared first in 
Deift [33], but we shall follow a somewhat different approach. We present 
the further results in supersymmetric quantum mechanics by recalling 
results obtained by B.Mielnik [27] concerning the search of superpartners 
for the harmonic oscillator. He pointed out that the factorization is not 
unique. He presented a new derivation of a important class of potentials 
previously obtained by P.B.Abraham and H.E.Moses with the 
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Gelfand-Levitan formalism [34]. Their energy and eigenfunctions can be 
directly obtained from the harmonic oscillator up to a zero mode state. In 
Section 3, we will show how this family of Hamiltonians is related to 
superintegrable systems with third order integrals. 

We consider the following Hamiltonian 

H„,r = \ . (2.9) 

We introduce the following first order operators 

a=^(-^ + x), a^ = ^{-^ + x) . (2.10) 
V2 dx v2 dx 

The Hamiltonian Hi and H2 are superpartners and have in fact the shape 
invariance properties 

a^a = Hosc-l = Hi, aa^ = Hose + I = ■ (2.11) 



This construction allows us to find the energy spectrum and the 
eigenfunction algebraically. B.Mielnik considered the Hamiltonian H2 and 
showed that the operators a and are not unique [27]. He defined the 
following new operators 



and required 



H2 = Hosc + l: = bb^ . (2.13) 



2 

He obtained the following Riccati equation [35] 

f3'{x) + p\x) = l+x'^ . (2.14) 

The fact of knowing a particular solution (/3(x) = x) allows to find the 
general solution [35]. He defined 

(3{x) = x + (p{x) , (2.15) 
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and found 



7 + Jq ax 

where 7 is a constant. There are two cases : with a singularity and without 
singularity. The inverted product h^h was not if2+const and was a new 
Hamiltonian 

H'^m^H..-,'i^)^--^.---i,^^j^^) . (2,17) 

We can obtain from H2 the creation and annihilation operators for H'. 
These operators are given by the following expression 

= b^a^b, s = b^ab , (2.18) 

with a and a"!" the annihilation and creation operators for H2. The 
eigenf unctions and energy spectrum of the Hamiltonian H' can be obtained 
from the Eq.(2.4). The coherent states have also been studied extensively 
[36]. This system is a special case of a one dimensional part of a 
Hamiltonian separable in Cartesian coordinates written in terms of the 
fourth Painleve transcendent. 



3 Higher order integrals of motion and 
SUSYQM 

Let us consider a two-dimensional Hamiltonian separable in Cartesian 
coordinates Ht{x, y, P^, Py) = Hx{x, P^) + Hy{y, Py) with creation and 
annihilation operators (polynomial in momenta) A^,, A\., Ay and A]^. These 
operators satisfy 

[//.,4] = A.4, [Hy,Al] = \yAl . (3.1) 
The following operators 

h = A^rAl f2 = A^Al- , (3.2) 
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commute with the Hamiltonian H 



[Ht, f^] = [Htj2] = 



(3.3) 



if 



mXx — nX. 



•y 



0, m,n E T, 



(3.4) 



Creation and annihilation operators allow us to construct polynomial 
integrals of motion. 

The following sums are also polynomial integrals that commute with the 
Hamiltonian H 



There are the integrals Ii and 12- The system Hf is thus superintegrable. 
By construction, the Hamiltonian Ht possesses a second order integral 
{K = Hx — Hy). The integral I2 is the commutator of Ji and K. The 
Hamiltonian Ht is thus superintegrable. We will show how supersymmetry 
makes it possible to construct superintegrable systems from 
one-dimensional Hamiltonian with creation and annihilation operators 
A^^ and A^. We choose in the y-axis a superpartner (or a family of 
superpartners) . This Hamiltonian Hy possess creation and annihilation 
operators that can be obtain from Eq.(2.7). A direct consequence of 
supersymmetry is the relation A^; = Aj^. We have thus the following integrals 



K = Hx-Hy, h = AlAy-AxAl h = AlAy + A^Al . (3.6) 



Let us, apply this construction to the interesting systems found by Mielnik. 
We take in the x axis the Hamiltonian H2 given by Eq.(2.9) and in the y 
axis its superpartner H' given by Eq.(2.17). We obtain a superintegrable 
system with integrals given by Eq.(3.6) with Eq.(2.10) and (2.18) 

K = H^- Hy, Ji = a\sy - ttxsl, I2 = a^sl + alsy , (3.7) 

where aj,, a^, sj^ and Sy are respectively the creation and annihilation 
operators of H2 and H'. 

These integrals are of order 2, 3 and 4. This superintegrable system appears 



T _ 4tm An _ Am a]""- T _ A\m An . Am a]i 
— ! -'2 — ^x ^y + ^x ^y 



n 



(3.5) 
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in the investigation of superintegrable systems with a second and a third 
order integrals separable in Cartesian coordinates. This is a particular case 
of a Hamiltonian written in terms of the fourth Painleve transcendent 
found by S. Gravel [23] and studied in Ref. 26. 



4 Construction of new superintegrable 
systems 

4.1 Hamiltonians involving the error function 

We consider the following superintegrable systems obtained in Ref. 23 and 
studied in Ref. 25 from the point of view of cubic algebras and SUSYQM 



We consider the case a = iao, G M. Let us define the two operators 



= + x^) + h'i^Z^ + TZ—X, + TZTT^]- (4-1) 



= ^{-h— + —^x + h{ — + — , 4.2 

V2 dx 2aQ X — ittQ x + la^ 

1 d ^ ^ ^ ^ 1 ^ "'^ )) (4 3) 

^/2 dx 2al X — iao x + iao 



We have 



2 ^ Sa^ + (x - zao)2 + (x + zao)^ + ia^ ' ^ ' 

The Hamiltonian Hg is the sum up to a constant of Hsi and Hs2- We apply 
the Mielnik's procedure to the Hamiltonian Hgi to find all the 
superpartners. We define the following operator 

^=^(^ + /3(^)), ^' = ^(-^+/5(^))' (4.6) 
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^0 = TTT^ + (:— :r + ) • (4-8) 



and demand Hgi = d)d. We obtain the following Riccati equation 

/3'(x)+/32(x) = |^ + — ^ + — ^ + ^ . (4.7) 

We know a particular solution 

1,1 1 

— + ( : \ ; 

Za^ X — ittQ X + zoo' 

Because we know a particular solution we can found the general solution. 
We consider 

/? = /3o(x) + , (4.9) 

and obtain the following the following equation 

(j)'{x) + (j)^{x) + 2/3o{x)(j){x) =0 . (4.10) 

We consider the transformation z{x) = and obtain a first order linear 
inhomogeneous equation 

- z\x) + 2/3o{x)z{x) + 1 = . (4.11) 

We obtain 

z{x) = e^o{al + x^f'y + -^{al + x^){2aox + e^v^(a^ + x^)Erf{-^)) 

(4.12) 

/3(a:) = -^x + ( — + — —)+ (4.13) 

zoq X — lao X + lao 

1 

e^(ag + ^2)27 + ^(a2 + x2)(2ao2; + e^o^{al + x'^)Erf{^)) 

Using the function z(x) given by the Eq.(4.12) the family of superpartner is 
thus given by 

i/, = i/,-^'W = | + g + ^-^ + ^-^ + g (4.14) 

_A[ \ , 

e^(a2 + x^Y-i + ^(ag + x^){2aox + e^v^(a2 + x^)Erf{^^) 
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The eigenf unctions and energy spectrum of Hamiltonian Hs^ have been 
obtained in Ref. 25 from supersymmetry. The eigenfunctions and energy 
spectrum of can be obtained directly from H^^ and Eq.(2.4). We can 
also obtain the creation and annihilation operators from those of Hgi. If we 
take Hx = Hgi and Hy = (the Hamiltonian is thus now given in 
term of the variable y) we obtain a new superintegrable Hamiltonian : 

^ ^^^^^ 

^ . 1 

dy 



+ 



^g'" \/2ao ' 



The creation and annihilation operators for the Hamiltonian Hgi are 

ml = clalcx, rrix = claxC^, (4.16) 

with 

The creation and annihilation operators of the Hamiltonian are 

rl = dlmldy, Ty = dlrUydy . (4.18) 

We can find from Eq.(3.6) the integrals of motion of the Hamiltonian of 
the order 2, 7 and 8 

K = Hx-Hy, Ii = miry - m^rl, I2 = mlvy + m^rl . (4.19) 

The integral I2 is given by the commutator of the integrals K and Ji. 
Because the harmonic oscillator is also isospectral to we have also the 
following superintegrable systems where we take H^. = and Hy = 

„ „ „ P2 p2 ^2^2 2 ^2 f^2 9^2 

Hf = Hx+Hy = + , + 4.20 

^ 2 2 8a^ Sal iv - ^^o) {y + laoY 4a^ 
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We have from Eq.(3.6) the following integrals of order 2, 5 and 6 

K = H^- Hy, Ii = airy - a^rl, h = a^Vy + a^rj. (4.21) 
The integral I2 is given by the commutator of the integrals K and Ji. 

4.2 Hamiltonians with fourth Painleve transcendent 

The following superintegrable system written in terms of the fourth 
Painleve transcendent can also be related to supersymmetric quantum 
mechanics [26] 

Hpi = ^ + ^+gi{x)+g2{y) , (4.22) 

9A^) = ^^Y^ ^\h^'^^ id j^x)+ujVnwxf{J -x) + — {-a+e) 

(4.23) 

92{y) = Yy" ■ (4-24) 

This Hamiltonian has a second and third order integrals. 

The function f is the fourth Painleve transcendent and /' = ^, z = ^JJ^x 

f{z) = P,{z,a,/3). (4.26) 



We will show that we can find new superintegrable systems from the 
superpartners of a one dimensional Hamiltonian with potential gi given by 
Eq.(4.23). This system was discussed in Ref. 26 and 30 and has a first and 
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second order super symmetry that allow to get the eigenfunctions and the 
energy spectrum. This system can have three, two or one infinite sequence 
of levels depending on parameters a and /3. When a potential possesses 
only one infinite sequence of energies, this potential may also allow singlet 
or doublet states. 



Let us consider 

H, = Pl^Y,{x) , i = \X (4.27) 
with a supersymmetry of order 1 and 2 with the following operators 

= ^{d^W{x)), q = -^{d^W{x)) , (4.28) 

= -2/i(x)<9 + 6(x), M = d'^ ^2h[x)d^h{x) . (4.29) 
From first order supersymmetry we have 

Vi = W'{x) + W\x), V2 = -W'{x) + W\x) - ^, (4.30) 

(another relations can be obtained from the supersymmetry of second 
order). The compatibility condition leads to 

W{x) = -h{x) - ^x . (4.31) 

The potentials Vi and V2 are obtained from (4.23) putting respectively 
e = -1 and e = 1 and adding hLj{^ - | - 1) ( with h{x) = y^fi^/jx)). 
We can apply the method to the Hamiltonian Hi{x) and find new operators 
k'' and k that factorize Hi. 

^ = ^(^ + /3(^)), ^' = ^(-^ + /5(^))' (4.32) 

This leads to a Riccati equation that we can solve because we know the 
particular solution W{x) and we find 

z{x) = -^ = e/^2^(-')'^^'(7+ re/^'^^(-")'^-"dx'), (4.33) 

0(X) J 

/3(x) = W{x) + . (4.34) 

Z[X) 
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We obtain 

Hsusy = ^ - ^im) (4.35) 

The eigenfunctions and the corresponding energy spectrum of Hi, H2 and 
thus Hp-^ were discussed in Ref. 26 and 30. Thus we can obtain directly with 
Eq.(2.4) eigenfunctions and corresponding energy spectrum of Hamiltonian 
Hsusy given by Eq.(4.35). We also know the creation and annihilation 
operators of the Hamiltonian Hi (and H2) and we can obtain from them 
the creation and annihilation operators for Hgusy by the supersymmetry. 
From these operators, we can form two integrals of motion and we have 
from the separation of variables in Cartesian coordinates an integral of 
order 2. This system is superintegrable. 



The creation and annihilation operators of Hi are given by the following 
third order operators 

at = q^M^ a = M^q , (4.36) 



with 

\f2h 



^i,2 = -oA/T/(A/T^)±(^^^^V^^^7^^if-^^) • (4-38) 



The creation and annihilation operators of Hgusy are given by 

t,t = k^a^k, V = k^ak . (4.39) 



The operators and v are quintic operators. If we take Hx{x) = Hi and 
Hy{y) = Hsusy we obtain the following Hamiltonian 

(4.40) 



13 



with the integrals of motion 

K = H^- Hy, Ii = a\vy - ttxvl, h = alvy + a^vl . (4.41) 

The integral I2 is given by the commutator of the integrals K and Ji. The 
integral K is of order 2, Ii is of order 7 and I2 of order 8. 



5 Conclusion 

In this article, we showed how supersymmetric quantum mechanics gives us 
a method of obtaining new superintegrable systems with higher order 
integrals of motion. Supersymmetry in quantum mechanics make it possible 
to find eigenfunctions and energy spectra from a superpartner using the Eq. 
(2.3) and (2.4). From a one-dimensional Hamiltonian and its superpartner 
we have constructed a two-dimensional superintegrable system and its 
integrals. The integrals are given by the Eq.(3.6). 

We discussed results obtained by B.Mielnik [27] in context of SUSYQM. We 
showed how we can generate a superintegrable system from the 
Hamiltonian he obtained and recover a particular case of a system with a 
third order integral from the Ref. 23 and studied in Ref. 26. 
From the method, we have explicitely constructed superintegrable systems 
written in terms of the error function and the fourth Painleve transcendent. 
These systems have higher integrals of motion. They possess respectively a 
second and a quintic integrals and a second and seventh order one. The 
supersymmetry allows also to find the wave functions and the energy 
spectrum. 

This method of generating superintegrable systems can be applied to other 
systems obtained in the contex of supersymmetric quantum mechanics. The 
results can be generalized in higher dimensions. 
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